Thermodynamics of a stochastic twin elevator 
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We study the non-equilibrium thermodynamics of a single particle with two available energy levels, 
in contact with a classical (Maxwell-Boltzmann) or quantum (Bose-Einstein) heat bath. The particle 
can undergo transitions between the levels via thermal activation or deactivation. The energy levels 
are alternately raised at a given rate regardless of occupation by the particle, maintaining a fixed 
energy gap equal to e between them. We explicitly calculate the work, heat and entropy production 
rates. The efficiency in both the classical and the quantum case goes to a limit between 100% and 
50% that depends on the relative rates of particle transitions and level elevation. In the classical 
problem we explicitly find the large deviation functions for heat, work, and internal energy. 

PACS numbers: 05.70.Ln,05.40.-a,05.20.-y 



I. INTRODUCTION 

Over the past decade, there has been growing inter- 
est in the stochastic energetics of small systems. This 
interest is driven in part by the impressive experimental 
and technological progress in bio- and nano-tcchnology. 
At the same time, the study of small scale systems has 
led to spectacular developments in nonequilibrium sta- 
tistical mechanics and thermodynamics. Brownian mo- 
tors and refrigerators [lj, work and fluctuation theo- 
rems 0, [1[ , and stochastic thermodynamics [3, Q provide 
prominent examples of these developments. Among the 
issues of specific interest arc the thermodynamic proper- 
ties of small scale stochastic systems and, in particular, 
the efficiency of interconverting different forms of energy. 
For classical heat engines, a certain degree of universal- 
ity has been identified for the transformation of heat into 
work. In particular, the efficiency rj at maximum power 
is found to be half of Carnot efficiency in the regime of 
linear response @ . This result has been illustrated by 
explicit calculations for several small scale engines 0- 
9]. However, many artificial and most biological engines 
operate in an isothermal environment. They transform 
one form of energy (e.g. chemical or electrical) into an- 
other form (e.g. mechanical or optical). Thermodynam- 
ics prescribes that the efficiency of this transformation 
is at most 100%, a limit again reached for a reversible, 
hence zero-power transformation. Concerning the effi- 
ciency at maximum power, it appears that there is again 
universality at the lowest order, i.e., in the regime of lin- 
ear response rj = 1/2 0. This result is reminiscent of 
the so-called maximum power transfer theorem from elec- 
trical engineering fill ] . enunciated by Moritz von Jacobi 
around 1840: maximum power is achieved when the load 
resistance is equal to the source resistance, with corre- 
sponding efficiency equal to 50%. The issue of universal- 
ity beyond linear response is currently under debate [lj| . 

The standard way to apply external work in statisti- 



cal mechanics is to systematically move (modulate) en- 
ergy levels or to modulate the potential energy. Unfortu- 
nately, the analytic treatment of even the simplest case, 
namely, modulating a single energy level, appears to be 
extremely difficult, see for example @, S El 0. The 
main purpose of this paper is to introduce an exactly 
solvable toy model which can be solved in full analytic 
detail, both in a classical and a quantum setting. We will 
call it the stochastic twin elevator. We present explicit re- 
sults for the rates of work, heat and entropy production. 
The efficiency of conversion of external work into internal 
energy is found to vary between 100% and 50%. For the 
case of a classical bath we also derive the analytic expres- 
sions for the large deviation functions [l5[ that charac- 
terize the statistics of the accumulated stochastic work, 
heat and internal energy and show that a heat fluctuation 
theorem is satisfied in the steady state. Such an explicit 
calculation is the exception rather than the rule [16j . 

The paper is organized as follows. In Sec. |IT]we intro- 
duce the stochastic twin elevator model and present the 
evolution equations and associated rates of heat, work, 
internal energy, and entropy production, as well as the 
results for the efficiency of the energy conversion process. 
In Sec. IIHI we concentrate on the classical bath and de- 
rive the steady state fluctuation theorem for heat. Wc 
also explicitly calculate the large deviation functions for 
the heat, the work, and the internal energy. We conclude 
with a brief summary in Sec. II VI 



II. STOCHASTIC TWIN ELEVATOR MODEL 

The model is defined as follows. A single particle, in 
contact with a heat bath, can reside in one of two avail- 
able energy levels separated by a fixed energy gap equal 
to e = e/kriT '. Here ks is the Boltzmann constant and T 
is the temperature of the bath. Thus e is the energy in 
units of /cbT; all energies will be expressed in these units. 
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The levels are alternately and instantaneously raised at 
random times at a rate k c , while maintaining the fixed 
energy gap e between them. When the particle occupies 
the level that is raised, the external agent must perform 
an amount of work equal to 2e on the system. If the 
level that is raised is unoccupied, the external work is 
zero. Note that this raising of the levels is a disturbance 
that drives the system away from equilibrium. Due to 
its contact with the heat bath, the particle can at any 
time make a thermal transition from the level it is occu- 
pying to the other level, absorbing from the bath (if the 
transition is uphill) or releasing to the bath (for downhill 
transitions) an amount of heat equal to s. In this way, 
work and heat can be monitored. At the same time, the 
entropy produced in the process is known from stochastic 
thermodynamics, see below for more details. 



The analysis of the properties of the system is thus re- 
duced to matrix algebra involving a time-independent 
4x4 matrix. 

Whether the bath is classical or quantum mechanical, 
the rates for thermal transitions of the particle obey the 
detailed balance condition 



(3) 



where W stands for the transition rate, h and I stand for 
higher and lower energy levels respectively. 



(a) 



The technical simplicity of the stochastic twin elevator 
is due to the following mapping onto a 4-state Markovian 
model. Let us arbitrarily call one of the elevator levels 
"1" and the other "2." The system can be in one of the 
four states {(1, 1), (1,2), (2, 1), (2, 2)}, the first index in- 
dicating which state is the lower and the second indicat- 
ing which state holds the particle. Thus, for example, 
the state label (1,2) means that level 1 is below level 2, 
and that the particle occupies level 2 (the higher one in 
this case). The system undergoes stochastic transitions 
between the four states. For example, lifting of the lower 
level in the state (1,1) corresponds to a transition to the 
state (2,1), cf. Fig. [T] Similarly, if the particle in state 
(1, 1) makes a thermal jump to the other energy level, the 
state changes to (1,2). These transitions all take place 
randomly in time, so that the probability distribution 
vector P with elements 



P(l,l)- 
P(l,2) 
P(2,l) 
P(2,2). 



evolves according to the master equation 
dV 



(1) 



dt 



= MP. 



(2) 



(b) 
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FIG. 1. Schematic representation of the two configurations of 
the twin elevator. The dotted line represents level 1 while the 
solid one is level 2. (a) Hopping of particle in configuration 1 
from the lower level to the higher level leading to a transition 
from state (1,1) to (1,2). (b) Transition from configuration 1 
to 2 due to lifting of the lower energy level by 2e thus changing 
the state from (1,1) to (2,1). 



A. Classical bath 

Consider first the case of a classical heat bath, and 
let k be the transition rate from the higher to the lower 
energy level in a given configuration. The transition rate 
from the lower to the higher is ke~ e (in general, k could 
depend on temperature): 



W(l,l)->-(l,2) 

W(i,2)-Ki,i) 



Wi 



(2,2) 



Wi 



(2,1) 



►(2,1) 
►(2,2) 



ke 
k. 



(4) 



As mentioned before, k c is the rate of lifting the lower 
energy level. We conclude that the transition matrix is 
given by the following expression: 



M 



-{ke- e + 
ke- £ 
k c 




-(k + k c ) 


kr. 



k c 
k c 
-(k c + k) ke~ e 

k -(ke- £ + k c ). 



= k 



-(e- E +0 1 £ 

e- £ -(1 + £ 

£ -(1 + e" e 

o e i -(e- e +e 



(5) 



where we have introduced the dimcnsionlcss parameter 
£ = k c /k. We will exhibit results in terms of the three 
parameters e, £, and k. We will focus on the steady 
state properties. The steady state probabilities P(l, 1), 



P(l, 2), P(2, 1), P(2, 2), with P(l, 1) + P(1, 2) + P(2, 1) + 
P(2,2) = 1, are found as the components of the right 
eigenvector of the matrix M associated with the zero 
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eigenvalue. Explicitly, 

P(1,1) = P(2,2) = 
P(1,2) = P(2,1) = 



e e (l+0 



2 [1 + e E (l + 2£)] ' 

l+ge e 
2[l + e £ (l + 2£)]' 



(6) 



In the absence of driving, k c = 0, the current between 
any two states vanishes. However, this is no longer so 
when k c > 0. In this case the four steady state currents 
between the states indicated by the subscripts are given 
by 



I — ^(1,1)-K1,2) — J(l,2)-K2,2) 

= P(l,l)far E -P(l,2)fc 

fcg(l -e £ ) 
" 2[l + e-(l + 20]' 



^(2,2)->-(2,l) — ^(2,1)-K1,1) 



(7) 



Let us look more closely at the nonequilibrium thermo- 
dynamic properties. The heat flux to the system (recall 
that all energies are given in units of fcgT) is given by 

Q= e[WW 2 P(l,l)-WWiP(l,2)] 
-e [W 21 ^ 22 P(2, 1) - IT 22 ^ 21 P(2, 2)] 
= ke [e^ (P(l,l) + P(2,2)) 

-(P(l,2)+P(2,l))]. (8) 

The rate of change of the work, that is, the power deliv- 
ered to the system, is 

W = 2e [W u _>ai-P(l, 1) + W 22 ^ 12 P(2, 2)] 

= 2efc£[P(l,l) + P(2,2)], (9) 

which reflects the fact that work in the amount of 2e is 
performed on the system when an occupied level is lifted. 
Using the First Law of thermodynamics, the increase in 
the internal energy of the system per unit time is just the 
sum of these two contributions, 



U = W + Q. 



(10) 



Finally, the rate of total entropy production associated 
with the master equation is given by [3] 



,P(v,j) log 



where the summation is over all possible states. 

The thermodynamic quantities Q, W, U and Si in non- 
dimensional form q, w, u and Sj can be rewritten as 



. = Q_ e(l - e E ) 
q ~ ki ~ l + e E (l + 2£)' 



w 



W _ 2£e E (l + Q 
k£ ~ l + e E (l + 2e) : 



(12) 



(13) 



U_ 
hi 



q = e, 



and 



(14) 



(15) 



We have plotted the results for the entropy production 
rate as a function of e and of £ in the two left panels 
of Fig. [2j Note that equilibrium can be reached in two 
different ways, namely, with e — > or with £ — > 0. In 
these limits s, goes to zero. When the system is out of 
equilibrium, Sj increases when e or £ increase. (The right 
panels will be discussed in the next subsection.) 




FIG. 2. (Color online) (a) Entropy production rate s; as a 
function of e for different values of £. The lines are for a 
classical bath and the symbols for a quantum bath, (b) s; 
versus £ for various values of e. The lindes are for a classical 
bath and the symbols for a quantum bath, (c) Difference in 
entropy production rates between the classical and quantum 
bath cases as a function of s for various values of f. (d) 
Difference in entropy production rates between the classical 
and quantum bath cases as a function of £ for various values 
of e. 



Finally, we turn to the efficiency of the system, which 
quantifies how efficiently the work done on the system is 
utilized in increasing its internal energy and is given as 



V 



U_ _ l + e £ (l + 2g) 
W ~ 2e e (l+£) 



(16) 



In the two left panels of Fig. [3] we have plotted the results 
for efficiency as a function of e and of £ (the right panels 
will be discussed in the next subsection). We observe 
that r\ decreases with increasing values of e, cf. Fig.^a), 
and increases with increasing values of £, cf. Fig. [3jb). 
We further note that: 

(i) For e — > 0, r\ — > 1. In this limit, all the four 
states have equal probability equal to 1/4, and so 
the efficiency r\ = 1. 
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FIG. 3. (Color online) (a) Efficiency r\ as a function of e for 
different values of £. The lines are for a classical bath and the 
symbols for a quantum bath, (b) r\ versus £ for different values 
of e. The lines are for a classical heat bath and the symbols 
for a quantum bath, (c) Difference of efficiencies between the 
classical and quantum bath cases as a function of e for various 
values of f. (d) Difference of efficiency between the classical 
and quantum bath cases as a function of £ for various values 
of e. 



(ii) In the other extreme limit, i.e., e — ¥ oo, rj 
1/2. In this limit 



'/ 



1 + 2g 
2(1 + 0' 



(17) 



which decreases from 1 when £ 
1/2 when £ 0. 



oo and goes to 



(iii) In general, there is a balance between the rate 
associated with configuration changes and that as- 
sociated with particle transitions. If the configura- 
tion changes very quickly compared to fc, the effi- 
ciency of the system increases. 



B. Quantum bath 

The difference between the classical and the quantum 
versions of our toy model lies in the nature of the bath. In 
the former, the bath is described by Maxwcll-Boltzmann 
statistics. In the latter, where, for example, the bath ex- 
citations might be phonons or photons, the statistics are 
Bose-Einstein. The rates at which the particle makes a 
transition between the two levels in a given configuration 
involve the emission or absorption of these excitations by 
the bath, and are now given by 



(i,i)-Ki,2) 



(2,2)->(2,l) 



kn(e) 



W(i,2)->-(i,i) = W / (2,i)->(2,2) =k[l + n(e)] , (18) 

where, as before, e is the energy difference between the 
levels in units of the thermal energy, n(e) = (e £ — 1) 1 is 
the Bosc-Einstcin distribution function, and k is a rate 
coefficient. Note that these transition elements obey the 
detailed balance condition Eq. ©. As in the classical 
case, the stochastic evolution of the system is described 
by the master equation Eq. ([2]), but now with the tran- 
sition matrix 



J 



M 



[n(e)+£] [n(s) + l] 

n(e) -[l + n(e)+€\ 
£ 

o e 





-[l + n(e)+£] 
He) + 1] 





n(e) 
Me) + 



(19) 



The master equation leads to the steady state solution 
for the probabilities 



P(1,1) = P(2,2) 
P(1,2) = P(2,1) 



e s (l + Q-£ 
2(2e e £ + e £ - 2£ + 1) ' 

e £ £ + 1 ~ j 
2(2e £ £ + e £ - 2£ + 1) ' 



(20) 



Following our earlier rules for the classical case, we can 
write the rates of heat, work, and internal energy influx 



into the system. In adimensional form we find 

e{l-e £ ) 
q ~ 2e £ £ + e £ - 2^ + 1 ? 
2 £ (e £ g + e E -Q 
2e £ £ + e £ - 2£ + 1 ' 
u = ill + q = e 

Si = -£q. 



w 



(21) 



In the left panels of Fig. [5] we show the entropy produc- 
tion rate as a function of e (upper panel, symbols) and 
of £ (lower panel, symbols). In the right hand panels of 
Fig. [2] we show the difference between the classical and 
quantum entropy production rates, in panel (c) as a func- 
tion of e for different values of £, and in panel (d) as a 
function of £ for different values of e. The classical and 
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quantum entropy production rates are equal in the lim- 
its e — > and e —> oo, but between these two limits the 
quantum entropy production rate is everywhere greater 
than in the classical case. As a function of £ for fixed e, 
the two again become equal as £ — > 0. As £ — > oo the 
difference goes to the limit — ee~ e /2. 

Finally, we calculate the efficiency of the system: 



U 2e e £ 



2£+ 1 



2(e^ 







(22) 



The efficiency as a function of e and of £ are shown by the 
symbols in the left hand panels of Fig. [3] We show the 
difference between the classical and quantum efficiencies 
(classical minus quantum) in the right hand panels of 
Fig. [3J for a number of parameter values. In the limits 
e —> and £ —> oo the difference between the two goes 
to zero, as seen in panel (c), as it should. In the limit 
£ — >• oo we find 



V = 



2£ + l 
2£ + 2' 



again exactly as in the classical case. The approach of 
the classical and quantum results to one another with 
increasing e is seen in panel (c) of the figure. For all 
values of £ between these limits the efficiency is higher 
in the classical case. As in the classical case, there is in 
the quantum case a balance between the rate associated 
with configuration changes and that associated with par- 
ticle transitions; if the configuration changes very quickly 
compared to k, the efficiency of the system increases. In 
panel (d) we show the difference between classical and 
quantum efficiencies as a function of £ for various values 
of e. Again, the efficiencies are equal when £ = 0, as they 
should be. As £ increases, the difference goes through a 
maximum (it is always higher in the classical case) and 
goes to zero again as £ — > oo and both efficiencies go to 
unity. 



III. FLUCTUATION THEOREM AND LARGE 
DEVIATION FUNCTIONS FOR CLASSICAL 
BATH 

We now turn to the classical case to explicitly calculate 
a number of other thermodynamic properties for the twin 
elevator system, in particular, the large deviation func- 
tion and the steady state fluctuation theorem for heat, as 
well as the large deviation functions for work and internal 
energy. 



A. Fluctuation theorem for Q 



Let Q be the accumulated heat transferred to the reser- 
voir up to time t. Since the transition of the particle be- 
tween the levels is stochastic, the total accumulated heat 
Q is also stochastic. Let P(i, Q, t) be the probability that 
the system is in state i at time t and the heat transferred 
to the reservoir is Q. The evolution of P(i,Q,t) follows 
from the equation 

P(i, Q,t + dt)=J2 Wj^idbPV, Q - AQj_w, t) 



I -^Wi^jdt 



P(i,Q,t). (24) 



Here AQj^ t is the amount of heat transferred to the 
reservoir as a result of a transition of the particle between 
the levels in state j, the resultant state being state i, with 
(23) rate Wj->i. As always, the quantities AQj-^i and Q are 
in units of fcgT. The differential form of the evolution 
equation follows from the master equation ([2]), 



j 

-J2Wi^P(i,Q,t). 



(25) 



We solve the above equation using the characteristic 
function (the subscript q labels this as the characteris- 
tic function for the heat), 



dQe~ XQ P(i,Q,t), 

-OO 



(26) 



whose evolution equation is obtained directly from 
Eq. (USD, 



dp q (i,X,t) 

at 



D -\AQj 



-^2Wi^jP q (i,\,t). 
3 

In matrix notation, it can be written as 

dp q 
dt 

where M g (A) is 



% = M g (A)p 9 , 



(27) 



(28) 



Mg(A) 



e^- 1 ^ -(1+0 

e o 

£ 







-Xe 





e (A-l)e 

-(e- E + 



(29) 
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and p q is a column matrix with four components. In 
writing the above matrix we make use of the fact that 
during an upward transition of the particle in a given 
configuration the reservoir loses heat to the system and 
so q is — £, and when the particle makes a downward 
transition heat flows to the reservoir, i.e., q = +e. We 
also note that the matrix M> and its adjoint satisfy 
the symmetry relation [l7j | 



M,(A) = Mj(l-A). 



(30) 



The matrix M 9 (A) has four eigenvalues, and the gen- 
eral solution of Eq. ([28)) is written as a linear combi- 
nation of the four independent associated eigenvectors. 
However, the large t behavior of p q is dominated by its 
largest eigenvalue, i.e., 



(31) 



where f q (X) is the largest eigenvalue of the matrix M 9 (A), 
which wc find to be 



ke e 



(-{2£e £ + e e + 1) + y^e 6 (f e e + e Ae + ed-^) + (1 + e £ ) 2 



(32) 



We observe that the eigenvalue f q (X) obeys the symmetry 
relation 



/,(A) = /,(!- A), 



(33) 



which is a direct consequence of the symmetry rela- 
tion ([3U| for the matrix M g (A). The above symmetry 
relation for the maximum eigenvalue reflects the steady 
state fluctuation theorem 0, i.e., 



Pig) 

P{-Q) 



(34) 



It is easy to verify that the average heat per unit time 
released to the reservoir is given in terms of the first 
derivative of the maximum eigenvalue evaluated at A = 0, 
i.e., 



(Q) 



df q (X ) 
dX 



A=0 



fc£e (e E 



1) 



l + e £ (l + 2£)' 



(35) 



We note that the magnitude of (Q) is the same as that 
given in Eq. (fT"2"|) (where the average is understood). 



We next turn to the explicit evaluation of the large 
deviation properties. For this we need to find the prob- 
ability P(Q,t) for long times. The heat Q is expected 
to grow linearly in time. We thus introduce the variable 
4> = Q/t = Q, which is the heat flux. The flux <fi can be 
positive as well as negative due to gain or loss of heat. 
According to large deviation theory [T{|, the probability 
P(Q,t) at long times can be written as 



P(Q,t)~e 



-tg q (<t>) 



(36) 



I 

where g q (4>) is the large deviation function. To find 
the relation between / g (A) and g q (4>), we implement the 
change of variables from Q to (f> in the integrand of 
Eq. I|3ip. which immediately allows us to identify f q (X) 
with the extrcmum of —(g q (<j)) + X4>) with respect to (f>. 
That is, invoking large deviation theory we see that the 
maximum eigenvalue f q (X) and the large deviation func- 
tion g q (4>) are related by a Legendre transformation, lead- 
ing to: 



3 9 (0) = -(/(A ) + 0A ), 
where is the solution of 

/^(A ) + 4> = 0, 



(37) 



(38) 



the prime denoting differentiation with respect to A. We 
find that the large deviation function g q (4>) is given by 
(see the Appendix) 



ke 



1) 



'4£ 2 e 2e + (1 + p/) 2 + 4£e £ 7 =F(0) + 



7 =F(0) 



--hi7 T W0, (39) 

£ 

where we have introduced the distinct notation g q (4>) = 
g q (4>) when (f> > and g q {4>) = fl^~(0) when <j> < 0, and 



7 T (0) = \ (p(4>) T v^(F^) , 



(40) 



with 



p{4>) 



A /4fc 4 ^ 2 e 4 e- £ + 4fc 2 e 2 2 ^ 2 + k 2 e 2 (j) 2 e-^(l 

k 2 e 2 £e~ E 



a e\2 



4</> 4 + 2</> 2 



(41) 
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In Fig. [4] we show the large deviation function and 
its derivative as a function of cj> for particular parameter 
choices. We note that the functions g£(4>) and gq{4>) 
match smoothly at <p = as their values are the same at 
that point, 



FIG. 4. (Color online) Large deviation function g q (<f>) and its 
derivative as a function of <f> f° r £ = 1, £ = 2, fc = 1. The 
solid line is for g^(4>) while the dashed line is for g~ ((/>). The 
value of (f) where g q (<f>) = is 0.23552, which is the same as 
(Q) obtained from Eq. ((35)) . 

I 



5 9 + (0) = 5 9 -(0) = -^fce- 



2£e £ 



1 - J8£e 3s / 2 + 4£ 2 e 2e + (1 + e e ) 



(42) 



as arc their first derivatives, 



dg+(0) _ dg~(0) _ 1 
dd> dd> 2 ' 



(43) 



We also note that g q (<p) has a single minimum at = 
(Q), as given by Eq. (|35[) . and at that point g q (4>) = 0. 
Although the large deviation function is a complicated 
nonlinear function of <p, the difference between g^(4>) and 
g q (4>) turns out to be a simple linear function. That is, 
for 4> > 0, 

g q (-4>) - 9M = g- q ■(-</>) - g+tf) = <t>, (44) 

which leads to the fluctuation theorem for heat as written 
in Eq. 



B. Large deviation function for VV 

We next turn to large deviation function for the work 
W done on the system. Work is done on the system 
whenever the filled lower energy level is lifted by energy 
2e. The process of lifting the level at rate k c is a random 
Markov process. As time increases, the work W also in- 
creases. Let P(i, W, t) be the probability that the system 
is in state i at time t with total accumulated work W. In 
analogy with the heat in Eq. (j2"5j) . the Markovian nature 

I 



I 

of the evolution allows us to write 
8P(i,W,t) 



dt 



Here AW 



= Wj^iP(j, W - AW H „ t) 

3 

-Y,Wi^P{i,W,t), (45) 

3 

2e is the work done on the system when 

z r . The 



its state change from j to i at rate W. 



J-H 



characteristic function is given by 

/>oo 

p w (i,X,t)= / dWe- xw P{i,W,t) 7 (46) 



with subscript w standing for work. The evolution equa- 
tion for p w (i, A, t) is 



d Pw {i,W,t) 
dt 



— AAW; 



i Pw(j, A, *) 



^ Wi-+jp w (i, X,t), 



or, equivalently, 



where 



dp w 
dt 



M w (X)p w , 



(47) 



(48) 



M w (X) = k 



-(e- £ +£) 1 £ 

+ o 

e" 2Ae £ -(1 + 

£ 1 



-2Ae, 



-(e- £ + 



(49) 
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Again, the large t behavior of p w is dominated by the where f w (X) is the largest eigenvalue of the matrix 
maximum eigenvalue, i.e., M^A), 



(e- xw ) 



(50) 



-ke 



l)e- A S/(4^ + l)e 2A£ + e 2 ( x +V s + 2e( 2A + 1 ) £ + 4e 2 <^ + 4£ 2 e 2£ 



(51) 



We can verify that the average work done per unit time, 
df w (0) 2fcee^(l + 



l + e £ (l + 2£)' 



(52) 



is the same as that obtained in Eq. ([13)) . 

Next we use the large deviation theory to shed more 
light on the probability distribution P(W,t), which al- 
lows one to write, 



I 

where g w {4>) is the large deviation function and </> is the 
flux of work (<fi = W/t = W). Here we note that (con- 
trary to the case of heat) <fi is always positive as the 
system always receives work so that W increases mono- 
tonically with time. The function gwifi) is again related 
to /to (A) by its Legendre transform. Following the same 
steps as for the heat, we now find the large deviation 
function gwifi) to be given by 



P(W,t) ~ e 



(53) 



9w{4>) = TTT 

2 d 



ke^[(2^ + e £ + l)Ci - Cz] - -Ci In 2d 2 

e 



(54) 



where 



Ci =<=' 



£(1 + -</> + yjk 2 t 2 e~ 2 ^ (4£e e + (1 + e £ ) 2 ) + 2 



(4^e e + (l + e e ) 2 ) 



(55) 



and 



A 



cb + ^k 2 e 2 e~ 2 ' (4£e e + (1 + e^) 2 ) + ^ 



(56) 



I 

In Fig. [5] we have plotted the large deviation function (subscript u for internal energy), as 
for work for a particular set of parameters. We see that 
it has a unique minimum and g w ((f) = at <j) — (W/t). 



dp u (i, U, t) 



iff 



\AUj 



(58) 



C. Large deviation function for the internal energy Here AU is the change in internal energy when the sys- 
tem undergoes a transition from one state to the other 

i . 1 1 ,ii j c c either due to a jump of the particle between the levels 

b mally, we turn to the large deviation function lor the . ,. r , . ,, n „ , , , , mi ^ , 

, „ ,, ,,° , or by lifting the tilled lower level, ihc above evolution 
internal energy following tnc procedure m the previous 

sections, we can write the evolution of the characteristic 

function, 



/oo 
dUe- xu P{i,U,t), 
-oo 



(57) 



equation can be rewritten as the matrix equation 

^=M u (X) Pu , 
where 



(59) 
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CD 




FIG. 5. (Color online) Large deviation function for work. 
Parameters are: e = 1, £ = 2, k = 1. The value of <f> where 
9w{4>) = is unique and its value is 0.23552, which is the 
same as that obtained from Eq. (|52|l . 



M„(A) = k 



e -(A+l)£ 
e 





-2Ae£ 



-(i + 
o 



-(i 





L + 

,Xz 







(60) 



In writing the above matrix, we have made use of the fact 
that during an upward transition of the particle, there is 
a gain of internal energy (AU = e) and, during a down- 
ward transition there is a loss (AU = —e). There are 
four such transitions where the system internal energy 
increases or decreases as a result of particle transitions 
between the levels. Similarly, when the lower level con- 
taining the particle is lifted, the internal energy of the 
system increases, AU — 2s. There are two such transi- 
tions that lead to an increase in the internal energy of 
the system by lifting. We are interested in the long time 
behavior, and so we can write 



-XU\ 



JMX) 



(61) 



where /«(A) is the largest eigenvalue of matrix M M (A) 
and is given as 



We verify that 



(62) 



(63) 



in agreement with Eq. (fl4")) . Next, we use large deviation 
theory to write for the probability distribution for the 
internal energy 



P{U,t) ~ e- t9 ^\ 



(64) 



I 

where g u (<l>) is the large deviation function for the in- 
ternal energy and cf> = U/t — U. The large deviation 
function and maximum eigenvalue are again related by a 
Legendre transform, from which we obtain 



9u{4>) = - ( 01og 



(65) 



In Fig. |6] we show the large deviation function for the 
internal energy for a particular set of parameter values. 
Again, we note that it has a unique minimum at 4> — (U) 
where g u {4>) = 0. 



IV. CONCLUSION 

We have presented the stochastic thermodynamics of 
a single particle that can reside on one of two energy lev- 
els, in contact with a classical or a quantum heat bath. 
The energy levels have a fixed energy separation and the 
system is driven out of equilibrium by alternately and 
stochastically lifting one of the two energy levels. The 
particle can make upward or downward transitions be- 
tween the energy levels mediated by the heat bath. At 
a given bath temperature, three parameters determine 
the behavior of the system: the energy separation be- 
tween the levels, the transition rates for the particle to 
move between levels, and the rate at which the levels are 
stochastically raised. The interest of this toy model lies 
in the fact that we can obtain explicit analytic expres- 
sions not only for average quantities such as work and 
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FIG. 6. (Color online) Large deviation function for the inter- 
nal energy. Parameters are: e = 1, £ = 2, k = 1. The value of 
4> where <7u(</>) = is unique and its value is 2, which is the 
same as that obtained from Eq. (|63[1 . 

heat flux, rate of entropy production and efficiency, but 
also of stochastic trajectory-dependent quantities includ- 
ing the steady state fluctuation theorem for the heat and 
the large deviation properties of work, heat and internal 
energy, both at the level of their characteristic function 
as well as in terms of the variables themselves. 



Appendix A: Large deviation function 

To find the large deviation function, using Eq. (|38|) we 
must solve the equation 



Differentiating f q (X) with respect to A gives 

df q {\) _ &k (e Xe -e' (A ' 1)£ ) 

(A2) 

Using Eq. (|A2[) in Eq. (|A1|) . we get the following 
quadratic equation in p: 

( 2 e 2 k 2 p 2 - ^e e <p 2 p - (4£ 2 £ 2 fcV + (f> 2 ^ 2 ) = 0, (A3) 

where 

p = l+— , 7 = e A6 . (A4) 

7 

p is always positive since it is a sum of exponentials. 
We thus consider only the positive solution of Eq. (|A3|) , 
which is given by Eq. (|41[) . We note that p is an even 
function of 4>, i.e., p(4>) = p(—<p), and is an increasing 
function of <f>. The limiting values of </> are, 

p{4>) = 2e e/2 for 0^0, 

= oo for <j) — > ±oo. (A5) 

Using Eq. (|A4[) . we get a quadratic equation in j((f>) 
with two roots 7 - (<^>) and 7 + (</>) as given by Eq. (|40[) . 
The valid solution for 7(0) must satisfy Eq. (|3"5)) . That 
is, for <^> > 0, df q (X)/d\ must be negative, which from 
Eq. (|A2[) requires that e^" 1 ^ < 1. Using the properties 
and limiting values of p(<fi) as given in Eq. (|A5|) . it turns 
out that this inequality is satisfied only if 7(0) = 7~(0). 
Similarly, for < 0, 7(0) = 7+ (</>). Finally, using these 
values of 7(0) in Eq. (f3"T)) . we get the desired large devi- 
ation function g q {4>) as given in Eq. (1391) . 
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